Photon bubbles in ultra-cold matter 
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We show that static and oscillating photon bubbles can be excited by diffused light in the laser 
cooled matter confined in a magneto-optical trap (MOT). The bubble instability is due to the 
coupling between the radiation field and the mean field oscillations of the ultra-cold gas, and it can 
provide a source for low frequency turbulence. We consider a diffusion dominated regime, which 
can be described by a radiation transport equation, coupled with the mean field equations for the 
cold atom gas. A perturbative analysis shows the occurrence of two different regimes with either 
oscillating or purely growing bubbles. This work could also be useful to understand similar processes 
in astrophysics. 
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Atomic physics has played an important role in the de- 
velopment of fundamental laws of quantum physics start- 
ing with the Bohr model for the atomic structure. With 
the advent of spectroscopy, laboratory based experiments 
in atomic physics have provided many important insights 
on astrophysical observations. As novel detection tech- 
niques in astrophysics allow for ever more subtle phenom- 
ena to be explored, ranging from quantum optics to las- 
ing in astrophysical systems [11, modern tools in atomic 
physics can be deployed to understand and predict pos- 
sible features in astrophysics. The recent study of pho- 
ton scattering in hot atomic vapors has thus lead to the 
direct observation of Levy flight of photons [2] induced 
by Doppler broadening and the possibility to implement 
dilute vapor based random lasing has been investigated 
[5]. Focusing on the impact of light on the atomic mo- 
tion, low frequency oscillations of the size of the MOT, 
reminiscent of oscillations found in variable stars called 
Cepheids, have been observed @] and numerical simu- 
lations predict turbulent like spatio-temporal structures 
[5]. The features are based on a repulsive force due to 
the exchange of scattered photons between nearby atoms 
[5J. This repulsive force can be described by an effective 
atomic charge, which explains the coulomb like expan- 
sion of the gas [7]. It can also explain the occurrence 
of collective processes of the plasma type, implying the 
existence of hybrid phonons, with a lower cut-off [HI [H], 
even though until now, these hybrid phonons have not 
been directly observed. 

Here we propose a new mechanism associated with the 
laser cooling process, which can lead to the formation of 
static and oscillating photon bubbles inside the gas, and 
eventually provide the source for the predicted turbulent- 
like spatio-temporal structures [5]. Photon bubbles have 
been considered in astrophysical context [TTJHT^] where 
huge photon densities are required to have any signifi- 
cant impact on high energy particles. The possibility of 
modeling astrophysical situations in laboratory experi- 
ments has been discussed in high energy physics [El [Mj . 
Due to the low kinetic energy of atoms accessible with 



the development of laser cooling techniques, radiation 
pressure effects can now be explored in laboratory based 
experiments with modest photon densities. Novel in- 
stabilities based on the coupling between atomic motion 
and light propagation come in the range of experimental 
observation. One class of such instabilities is based on 
"atom-field" coupling, where the index of refraction of 
the atomic density affects the propagation of light which 
couples back to atomic motion [T5lll8| . This regime can 
be explored in clouds of cold atoms with low optical thick- 
ness, using e.g. large atom-laser detuning. A second 
regime is based on "atom-intensity" coupling, where co- 
herent effects and interferences can be neglected. This 
regime valid for larger optical thickness will be considered 
in this letter. Note that in contrast to past observation 
of pattern formation in hot atomic vapors [THl HO] the 
instabilities based on the coupling to the atomic motion 
do not require non linear optical response of the atoms. 

Our model is based on the coupled photon and atom 
density evolution equations. The standard approach in 
the laser cooling community is to distinguish two terms. 
The first term is due to the coherent coupling between 
the incident laser beams and the atoms. This term de- 
scribes the friction and trapping terms producing the 
confined cold atom sample. A second term, aimed at 
describing the collective, many atom case, is based on 
the attenuation of the incident laser beams, often called 
the shadow effect [3T] , and of multiple scattering between 
the atoms [5J . In MOT experiments with very large num- 
ber of atoms, the multiple scattering (diffusion) terms 
dominates over the shadow effect. The combined atom- 
light interaction contain terms proportional to the atom 
number density. Therefore, a perturbation in the density 
will lead to a perturbation in the local photon intensity, 
which can become unstable. If the instability is uniform 
and isotropic, a local bubble can be formed, where both 
the photon intensity and atom density are modified. 

Let us consider the most relevant situation for the 
MOT conditions, where the incident laser cooling beams 
are deflected by multiple scattering, randomizing the 
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photon propagation direction, and the diffusion effects 
become dominant. In order to describe this new situa- 
tion we start from the energy transport equation for the 
photon field which, in a region where radiation sources 
are absent, can be generally stated as 

d 

— I u + V-S u = -i u I u . (1) 
at 

Here, 7^ describes photon absorption at the relevant 
wavelength by e.g. hyperfine Raman scattering. Even 
though the loss can often be neglected in realistic de- 
scription of experimental situations, we keep it in our 
model as a finite non zero loss term can induce a finite 
threshold for the instabilities predicted in this letter. The 
quantity I u can be defined as the spectral energy density 
of the electromagnetic radiation W(u>, k), integrated over 
all the possible directions of propagation. If radiation is 
made isotropic by multiple scattering, the energy flux is 
determined by a diffusive process, characterized by (see 
e.g. [23 |23]) S w = -DVI U and lu] = Dk 2 a , where k a 
is the inverse of the energy absorption length. The dif- 
fusion coefficient is determined by D — 1 2 /t, where the 
photon mean free path is I — l/n a aL, where n a is the 
atom number density, and the laser atom scattering 
cross section. The photon diffusion time r can be con- 
sidered as nearly independent from the atom density, as 
shown by cold atom experiments [24]. We therefore get 
D oc n~ 2 . 

On the other hand, the atom density n a can be deter- 
mined in the mean field approximation by the following 
fluid equations, which determine the atom density n a and 
mean velocity v, 

dn a „ , s n dv , ^ N F VP 

(2) 

where P is the pressure, v is the damping rate result- 
ing from the viscosity of the gas, and M is the atom 
mass. The collective force F, resulting from the ex- 
change of photons between nearby atoms is determined 
by a Poisson type of equation, V • F = Qn a , where 
Q = ((tr.-(Tl)o'lIuj defines an effective charge for a single 
atom in the mean field, and an and <jl are the atom radi- 
ation scattering and atom laser absorption cross sections 
[11 [7]. We notice that the charge parameter Q is pro- 
portional to the laser intensity, thus providing a coupling 
between the atom density n a and the photon intensity I u . 
In contrast with the astrophysical models for the dynam- 
ics of the gas around neutron stars |10j , young high-mass 
stars [12] or black hole accretion disks [11] , the fluid equa- 
tions for the atomic distribution in our situation do not 
depend on magnetic fields, preventing a direct mapping 
of the predictions for photon bubbles in astrophysics. 

The two coupled eqs. (1) and (2) provide the starting 
point for our analysis. Using a perturbation analysis, 
we assume that I u = Iq + I, n a = tiq + h, and treat v 



and F as perturbations. Linearizing these equations with 
respect to the perturbed quantities, and assuming that 
they evolve asexp(iq-r — Hit), we arrive at the dispersion 
relation 

(-ifi + D q 2 + 7o ) (fl 2 + iSlv - to 2 ) = -j8(e + ia)Q, (3) 

where we have used e = gV 2 Jo, a = g(q • VJo), and 
uj 2 = uj 2 +u 2 q 2 , with the effective plasma frequency uj p = 
(Qno/M) 1 / 2 , and the sound speed u s [5]- We have also 
used the unperturbed diffusion coefficient Dq — \/n 2 a 2 L T, 
70 = D k 2 , and g = 2Dq/tiq. Finally, we have introduced 
the coupling coefficient (3 = w 2 no/io- 

At this point it should be noticed that, in the ab- 
sence of coupling between the photons and the atomic gas 
(/3 = 0), this dispersion relation would describe two in- 
dependent modes, a purely decaying photon mode, such 
that f2 = —i(D q 2 + 70), where damping results from 
both diffusion and absorption, and an oscillating fluid 
mode, determined by f2 = uj s — iis, which decays due 
to viscosity. The dispersive properties of this oscillating 
mode have already been discussed in detail, for both the 
fluid and kinetic regimes [5] . It is now interesting to con- 
sider the modes which result from the coupling between 
the photons and the atomic g clS, (IS described by the new 
dispersion equation ([3]), when j3 7^ 0. 

Let us first examine oscillating perturbations with fre- 
quency of order ui s . Using fl = lo s + S, and assuming 
\5\ lo s , we obtain solutions of the form 

6=-L-M^(i + °^ + *). (4) 
2 \ lo s ui s J 

This can be considered as the perturbed oscillating fluid 
mode. We see that, apart from a small correction in 
the real part of the frequency, we can have a growth or 
damping rate, determined by the quantity T = 3(5) = 
Sy(f2). We focus on the case where the quantity a can be 
neglected and the following inequality is verified 

|q-V7 |« ^g 2 V 2 / . (5) 

In this case, instability will occur if the following two 
conditions are simultaneously satisfied 

/?e < , > v. (6) 

The first condition implies that, for j3 > (i.e. for 
> &l) oscillating bubbles can only occur in regions 
with V 2 io < 0. Notice that a positive /3 corresponds to 
the most natural experimental scenario in a MOT. The 
sign of V 2 Jo depends on boundary conditions or fluctu- 
ations of the unperturbed light intensity (e.g. due to 
imperfections) and can vary across the cloud. The emer- 
gence of instabilities thus depends on local initial condi- 
tions. We note that in configurations where the shadow 
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FIG. 1: Schematic representation of spherical bubbles result- 
ing from the unstable coupling of laser light with the atom 
mean Geld. Intensity plots of the function F(r), as determined 
by eqs. pi). 



effect would dominate over the repulsion, with j3 < (i.e. 
&r < instabilities can occur in regions of V 2 /o > 0. 
On the other hand, the second condition in Q deter- 
mines the instability threshold, where the growth rate 
has to compensate for the losses due to the viscosity in 
the ultra-cold gas. In the absence of viscosity, or for 
conditions well above the threshold, the growth rate at- 
tains its maximum value T max = \/3e\/2u> 2 . This will lead 
to the formation of photon bubbles, as discussed below. 
In the opposite situation where a becomes the dominant 
term in Q and the inequality |5]) is reversed, we still 
have mode instability if q • VIq < 0, but isotropy is lost 
and the growth rate will maximize along the gradient of 
the photon intensity. 

Let us now turn to purely damped (or growing) pertur- 
bations, with 5R(il) ~ 0. This can now be considered as 
the perturbed photon decay mode. Assuming now that 
|£!| 2 <C w 2 , equation ^ leads to 



r = 



^ - (^g 2 + 7o ) 



(7) 



This shows that unstable modes can indeed exist, for 
(3e > 0. It should also be noticed that these modes are not 
exactly purely growing modes, because a residual value of 
the mode frequency still exists, 5i(Sl) ~ (aD Q /2I ). This 
quantity goes to zero with a = q- Vio, and therefore van- 
ishes in the spherically symmetric case. For commodity, 
we will keep calling these modes, purely growing modes. 

It is useful to compare these modes with the above 
oscillating modes. First, they occur for opposite signs of 
the quantity /3e, the oscillating modes for a negative sign, 
and the purely growing modes for a positive sign. The 
maximum growth rate for the purely growing modes are 
two times larger than for the oscillating modes, T max = 



fie/uj 2 . On the other hand, the threshold conditions are 
different: i) For the oscillating modes, eq. ^ shows that 
the instability is limited by the viscosity of the atomic 
gas. ii) In contrast, for purely growing modes, eq. Q 
shows that the instability is limited by photon diffusion 
and losses. 

In this comparison we have assumed a negligible a, 
when the instabilities can be seen as isotropic. In this 
situation, we can associate such unstable modes with 
the formation of photon bubbles. In order to describe 
more explicitly the formation of such bubbles, we replace 
the plane wave modes by spherically symmetric perturba- 
tions described by (1, fi) on F(r) exp(— iQt), where F (r) 
satisfies the equation 



V 2 F(r) = 



r 2 dr 



c) 



F(r) = ~q 2 F(r), (8) 



where q 2 is a positive real quantity. Going back to 
the perturbed equations, but ignoring the quantity a 
q • Vio which would brake the spherical symmetry, we ar- 
rive again at the dispersion relation ([3]), but with a = 0. 
Therefore, our previous analysis for the oscillating and 
purely growing modes remains valid for such spherically 
symmetric perturbations. Equation ^ can be written 
as a spherical Bessel function, with a well known non- 
singular solution, given by 



F(r) 



Ji/2{qr) 



sin(qr) 



qr 



(9) 



where Ji/2(qr) is a Bessel function of the first kind. A 
linear superposition of four such solutions with arbitrary 
amplitudes and random positions is represented in Fig. 
[T] for illustration. Of course, superposition brakes down 
in the nonlinear regime, where instability saturation will 
occur. For purely growing modes, the saturation level can 
be easily estimated by replacing in the threshold condi- 
tion the equilibrium diffusion coefficient D by its 
perturbed value Z?o n o/( n o ~ l r ~ l l) 2 - Saturation will then 
occur for 



Msat — n 



D Q q 2 



7o 



1/2 



(10) 



Saturation can then be understood as a result of an in- 
crease of the diffusion losses due to a local density de- 
pletion. According to the coupled equations, a local de- 
crease in the atomic number density h also leads to an 
increase of the local photon intensity. Near the instabil- 
ity saturation, 3(f2) ~ 0, for the static bubbles we have 
the following relation I — —n/(Doq 2 + 70). Therefore, 
we have opposite signs for the perturbations I and n, as 
expected. 

Finally, it is useful to write the threshold conditions 
and growth rates in a more explicit form, in terms of the 
typical time and scale lengths. Introducing the photon 
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FIG. 2: Threshold curves for photon bubble formation in the 
ultra-cold gas: (a) Oscillating bubbles, photon characteristic 
scale L/L n , versus normalized wavenumber z — q/k a with 
L n — (uj p /uj B )y/2/k a ; (b) Purely growing bubbles, L/L n ver- 
sus z = \frv, with normalization factor L n = (u! p /lo s )IV^- 
Instability occurs below theses curves. 



intensity scale length L, such that L~ 2 = Iq \V 2 Iq\, we 
can write the threshold and growth rate for the oscillating 
bubbles as 



L < 



■21 



1P_ 
tL 2 ' 



(11) 



A similar analysis can be made for purely growing bub- 
bles, leading to 



L < 



0J S (q 2 + fc2)l/2 



2P_ 
tL 2 ' 



(12) 



Given that uj p ~ uj s , and assuming k 2 3> q 2 , we get a 
threshold condition L < l/k a , which can be easily satis- 
fied. On the other hand, taking the experimental condi- 
tions of [23], we have / ~ 300/im, D ~ 0.66m 2 /s, corre- 
sponding to t ~ 0.1/xs. Thus, even for large L ^> I, the 
growth rates T max can stay well above ls^ 1 . This can be 
satisfied in current experimental conditions, with oppo- 
site signs of /3e, for the oscillating and the purely growing 
cases. The threshold conditions are represented in Fig. 
2. Instability occurs for an intensity scale length L ly- 
ing below the curve (a) for oscillating bubbles, and below 
curve (b) for purely growing ones, as determined by eqs. 
(Ill and (12) respectively. For large wavelengths, such 



that qL < 1, the above results are only approximate, and 
our local analysis should be replaced by a global mode 
analysis. 

In conclusion, we have studied the stability of low fre- 
quency perturbations in the photon intensity and atom 
number density in the laser cooled gas. We have con- 
sidered the diffusive regime, where multiple scattering 
of photons by the cold atoms randomizes the direction 



of photon propagation. Similar to astrophysical systems, 
radiation pressure is also expected to be important to cre- 
ate turbulence in cold atomic vapors by destabilizing the 
atomic density distribution with subsequent feedback on 
radiation transport. We have used the photon transport 
equation, coupled to the fluid equations for the ultra-cold 
gas. Such a coupling is mediated by the atom effective 
charge, which depends on the photon intensity, and by 
the diffusion coefficient, which is inversely proportional 
to the square of atom number density. 

We have shown that, both static (or purely growing) 
and dynamic (or oscillating) photon bubbles can be ex- 
cited. These two distinct modes result from a pertur- 
bation of the two natural modes of the two fluids, the 
photon gas and the atomic gas, when considered sepa- 
rately. Due to coupling between the two interacting flu- 
ids, mediated by both the diffusion coefficient and by the 
effective charge, the two natural modes eventually be- 
come unstable, leading to the formation of purely grow- 
ing and oscillating structures. In certain experimental 
conditions, such that the quantity a = q • V/ ~ 0, these 
unstable structures acquire a spherical symmetry. In con- 
trast, when 0, new oscillating modes appear, which 
depend on the direction of propagation and which grow 
faster along the gradient of the photon intensity. We 
have therefore characterized the conditions under which 
photon bubbles can be excited. Nonlinear saturation lev- 
els have also been estimated. As shown above, photon 
bubbles represent a local increase of the photon num- 
ber density, associated with a local decrease of the atom 
number density. 

The static and the oscillating bubbles occur for dif- 
ferent signs of the quantity (3e, which is typically de- 
termined by the sign of V 2 2o- Such bubble instabilities 
are isotropic, and a simple spherically symmetric solu- 
tion was identified, although the basic ingredients of this 
intensity-atom coupling induced instabilities should also 
occur in different geometries. These instabilities could 
provide the driving mechanism for random structures 
and the excitation of low frequency turbulence in cur- 
rent MOT experiments. Experimental investigations of 
these structures are in reach with existing fast imaging 
techniques and we expect various geometries and excita- 
tion techniques to provide experimental insight for these 
instabilities. Finally, it should be added that radiation 
pressure is also expected to be at the origin of turbulence 
in astrophysics, opening the way to simulate complex as- 
trophysical situations in the laboratory. 
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